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The  objective  of  this  program  is  to  conduct  a  study,  design  and  fabricate 
a  cold  crucible  system  for  the  synthesis  and  Czochraiski  crystal  growth  of 
refractory  semiconductors  and  oxides.  During  the  first  6  months  a  theo¬ 
retical  analysis  of  the  cold  crucible  relating  the  proper' ies  of  the 
material  to  the  electrical  characteristics  if  the  cold  crucible  has  been 
essentially  completed.  The  thermal  analysis,  which  includes  thermal 
modeling  by  means  of  a  computer  technique,  has  also  been  completed.  Ar 
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initial  detailed  analysis  of  the  RF  power  control  system  necessary  to 
provide  the  desired  stability  for  crystal  growth  has  been  completed. 

The  design  of  a  cold  crucible  assembly  suitable  for  the  growth  or  re¬ 
fractory  oxide  crystals  has  been  completed.  It  is  a  simple,  rugged, 
j  modular  assembly  which  will  be  readily  adaptable  in  the  future  to  a 
control led-atmosphere  furnace  system. 


We  plan  to  use  a  50  KW  (5-?  MHz)  PF  power  generator  manufacfured  by 
Fritz  Heuttinger  Elektronik  GmbH  (Freiburg,  Germany)  for  the  performance 
testing  of  the  cold  crucible  assembly.  This  particular  system  was 
selected  because  of  its  unique  circuit  design  to  provide  for  high 
stability  under  widely-varying  load  conditions.  Delivery  of  this  system 
from  Germany  is  expected  shortly. 


TECHNICAL  SUMMARY 


The  objective  of  this  program  Is  to  conduct  a  study,  design  and  fabricate 
a  cold  crucible  system  for  the  synthesis  and  Czochralskl  crystal  growth  of 
refractory  semiconductors  and  oxides.  Specifically,  the  work  to  be  per¬ 
formed  is: 

Study  by  theoretical  analysis  of  the  cold  crucible  to  relate  the 
properties  of  materials  to  the  electrical  and  thermal  character! sties 
of  the  cold  crucible  and  define  a  control  system  to  provide  a  high 
degree  of  stability  for  crystal  growth. 

Design  a  cold  crucible  assembly  complete  with  control  system  suitable 
for  oxide  crystal  growing.  Define  the  criteria  necessary  for  inte¬ 
grating  the  cold  crucible  system  with  a  controlled  atmosphere  (up  to 
100  atmospheres)  crystal  growing  furnace. 

Fabricate  a  cold  crucible  assembly  complete  with  control  system. 

Test  and  evaluate  these  to  explore  the  operating  techniques  suitable 
for  crystal  growing  applications,  and  demonstrate  high  temperature 
operation  under  oxidizing  conditions. 

During  +he  past  6  months  the  following  topics  were  covered: 

I.  The  theoretical  analysis  of  the  cold  crucible  relating  the  properties 
of  the  material  to  the  electrical  characteristics  of  the  cold  crucible 
was  completed.  The  thermal  analysis,  which  includes  thermal  mooeling 
by  means  of  a  computer  technique,  has  also  been  completed.  During 
the  course  of  the  program  it  Is  anticipated  that  it  may  prove 
necessary  to  refine  the  analyses  following  the  operating  tests  of 
the  cold  crucible  assembly. 
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2. 


A  detailed  analysis  of  the  RF  power  control  system  necessary  tc  provide 
the  desired  stability  for  crystal  growth  has  been  completed.  It  may 
also  be  necessary  to  partially  modify  the  control  system  following  runs 
carried  out  in  the  cold  crucible  assembly. 

3.  The  design  of  a  cold  crucible  assembly  suitable  for  the  growth  of  re¬ 
fractory  oxide  crystals  has  been  completed.  Based  upon  the  guidelines 
established  by  the  theoretical  analysis,  as  well  as  first-hand  know¬ 
ledge  of  cold  crucible  configurations  which  are  currently  operating 
with  considerable  success  in  +he  U.S.S.R.,  we  have  designed  a  simple, 
rugged,  modular  assembly  which  will  be  readily  adaptable  in  the  future 
to  a  control ied-atmosphere  furnace  system. 

4.  We  plan  to  use  a  50  KW  (5-7  MHz)  RF  power  generator  manufactured  by 
Fritz  Huettinger  Elektronik  GmbH  (Freiburg,  Germany)  for  the  perfor¬ 
mance  testing  of  the  cold  crucible  assembly.  This  particular  system 
was  selected  because  of  its  unique  circuit  design  which  provides  for 
high  stability  under  widely-varying  load  conditions. 

Celivery  of  this  system  from  Germany  was  originally  expected  by  year 
end  1973;  however,  iT  has  been  seriously  delayed  due  to  t:.9  world-wide 
shortage  of  vendor-supplied  electronic  components.  We  have  recently 
been  informed  that  the  apparatus  will  be  undergoing  testing  during  the 
last  week  of  April  and  hopeful iy  will  be  air-shipped  to  the  United 
States  by  the  end  of  May. 
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PREFACE 


Major  advances  in  the  development  of  cold  crucibles,  i.e.  vater-cooled  metal 
crucibles  for  use  with  RF  induction  heatirg,  used  for  "skull  melting"  have 
taken  place  in  the  past  two  years  with  significant  impact  on  refractory 
materials  technology.  For  the  first  time,  it  is  possible  to  contain  stable, 
un contaminated  melts  of  even  the  most  refractory  materials  tor  which  no 
crucible  exists. 

Russian  scientists  at  the  Lebedev  Physical  Institute  in  Moscow  have  developed 
skull  nelting  to  a  point  where  it  has  become  a  reproducible  production 
opera* ion.  Oxide  melts  at  temperatures  up  to  3G00°C  with  weights  up  to 
25  ki icgrams  are  being  used  in  the  production  of  a  wide  variety  of  strategic 
materials  including  unique  laser  crystals,  refractory  optical  elements  and 
melt-cast  ceramics. 

The  new  process  involves  direct  high-frequency  induction  heating  of  the 
material  contained  in  a  water-cooled,  crucible-like  structure.  The  melt 
fooned  is  contained  by  a  sintered  shell  or  "skull"  of  identical  composition 
so  that  the  problems  of  reaction  and  contamination,  traditionally  the  :nost 
severe  problem  faced  in  the  containment  of  refractory  melts,  have  been 
virtually  Jliminated.  The  melt  can  then  be  cast  or  recrystallized  using 
standard  methods  of  crystal  growth. 

During  the  past  few  years,  high-pressure  crystal  growing  furnace  systems 

0 

developed  in  the  United  Stages  by  Arthur  D.  Little,  Inc.,  have  enabled  the 
materials  scientist  to  cope  with  the  problems  of  high  vapor  pressure  in 
producing  single  crystals  of  new  quantum-electronics  materials  such  as 
Gallium  PhospLide.  It  now  appears  feasible  to  combine  the  cola  crucible 
with  its  inert  ultra-high  temperature  capabilities  with  domestic  high- 
pressure  technology  to  provide  a  maferials  research  and  production  tools 
of  unprecedented  capabilities. 
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This  report  stalls  work  performed  in  the  first  6  months  of  the  program 
during  which  time  theoretical  studies  were  completed,  design  of  the  cold 
crucible  assembly  finished  and  work  initiated  on  Its  construction. 
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1.0  INTRODUCTION 


The  principal  objective  of  the  proposed  program  is  the  design  and  develop¬ 
ment  of  a  skull  melting  system  which  can  be  used  fcr  the  production  of  high 
purity  single  crystal  oxides.  Typically,  most  of  the  cold  crucibles  in 
operation  today  have  evolved  by  empirical  trial  and  error  procedures.  A 
theoretical  analysis  of  the  cold  crucible  has  been  carried  out  to  relate 
the  properties  of  the  material(s)  to  be  melted  to  the  electrical  and  thermal 
characteristics  of  the  cold  crucible.  Based  upon  tnis  analysis  we  have 
established  basic  guidelines  for  the  design  and  construction  of  cold  crucible 
assembl ies. 

One  of  the  major  problems  as  yet  unresolved  in  the  use  of  the  cold  crucible 
for  the  growth  of  single  crystals  is  related  to  the  stabilization  and  pre¬ 
cise  control  of  melt  temperature.  We  have  analyzed  the  critical  control 
parameters  of  the  cold  crucible  and,  based  upon  the  unique  RF  power  control 
techniques  which  have  been  developed  in  the  United  States,  have  defined  a 
control  system  to  provide  the  high  degree  of  stability  necessary  for  crystal 
growing  applications. 

A  cold  crucible  assembly  complete  with  control  system  suitable  for  oxide 
crystal  growing  operation  has  been  designed  and  is  being  constructed. 

Testing  and  evaluation  of  the  cold  crucible  assembly  and  control  system  will 
be  carried  out  to  explore  the  operating  techniques  suitable  for  crystal 
growing  applications.  An  oxide  such  as  stabilized  zirconta  (M.P.  2690°C) 
will  be  used  to  demonstrate  very  high  temperature  operation  under  oxidizing 
conditions.  We  will  then  deliver  a  cold  crucible  assembly  to  the  Air  Force 
Cambridge  Research  Laboratories. 

The  ultimate  goal  of  this  work  is  the  integration  of  the  cold  crucible  system 
into  a  crystal  growing  furnace  chamber  capable  of  operation  under  control  I 'd 
atmosphere  at  pressures  up  to  100  atmos.  It  should  be  realized,  however. 


8 


2.0  TECHNICAL  DISCUSSION 


2.1  Introduction 


During  the  past  6  months  we  have: 

1.  Theoretically  analyzed  the  electrical  characteristics  of  the 
cold  crucible  system.  The  problems  involved  and  the  conclusions 
are  presented  in  Section  2.2  while  the  derivation  is  discussed 
in  Appendix  A. 

2.  Theoretically  analyzed  the  thermal  characteristics  of  the  cold 
crucible.  The  conclusions  are  presented  in  Section  2.3  and  the 
derivation  shown  in  Appendix  B. 

3.  Theoretically  analyzed  the  control  of  the  RF  power.  This  is 
discussed  in  Section  2.4  and  the  derivation  is  presented  in 
Appendix  C. 

4.  Designed  and  initiated  construction  on  the  cold  crucible  assembly. 
Refer  to  Section  2.5  for  review  of  the  design  criteria. 

2.2  Theoretical  Analysis  of  the  Electrical  Characteristics  of  the  Cold 

Crucible 


The  following  problems  involved  in  skull  melting  of  refractory 
materials  have  been  theoretically  analyzed: 

•  The  effect  of  the  copper  cage  on  the  magnetic  field  applied  to 
the  refractory  charge. 

•  The  resistive  and  Inductive  effect  of  the  cage  on  the  radio 
frequency  coil. 
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°  The  effect  of  the  skin  depth  of  field  penetration  on  the  total 
power  input  to  the  charge. 

°  The  spa+!ii  Jlct'-ibution  of  the  power  input. 

tf  The  resistive  loading  of  the  RF  coil  due  to  power  absorbed  by 
the  charge. 

#  The  RF  magnetic  field  required  to  melt  metal  particles  in  the 
charge  to  achieve  startup. 

it  is  found  that  a  properly  designed  cage  has  no  effect  on  the  field 
applied  to  the  charge.  The  cage  merely  lowers  the  inductance  of  the 
RF  coil  slightly  by  excluding  the  field  from  the  volume  occupied  by 
the  cage  tubing  and  adds  a  small  resistive  loading  to  the  coil. 

The  power  delivered  to  the  molten  charge  by  a  given  applied  RF  field 
is  calculated  for  a  disk-shaped  charge  approximated  by  an  ellipsoid 
of  revolution.  This  is  done  for  the  field  penetration  skin  depth  both 
large  and  small  compared  with  the  smallest  semi-axis  of  the  ellipsoid. 
These  two  limiting  cases  converge  toward  a  common  maximum  value  of  the 
power  when  the  skin  depth  in  each  case  approaches  the  size  of  the 
semi-axis.  For  an  ellipsoid  of  dimensions  8cmx8cmx2cm  in  a 
field  of  1000  A/m,  this  maximum  power  is  several  kilowatts. 

The  calculations  show  that  the  heating  power  is  strongly  concentrated 
around  the  equator  of  the  ellipsoid,  especially  when  the  skin  depth 
is  small  compared  with  the  thickness  of  the  ellipsoid.  An  expression 
is  derived  for  the  resistive  loading  of  the  RF  coil  by  the  induced 
currents  in  the  molten  charge. 
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Startup  by  means  of  2  mm  diameter  metal  spheres  embedded  in  the 
powdered  charge  is  found  to  be  theoretically  feasible  but  requires 
a  starting  field  perhaps  10  times  larger  than  that  needed  to  maintain 
the  melt. 

2.3  Thermal  Analysis  of  the  Cold  Crucible 

During  this  quarter  the  thermcl  analysis  of  the  melt  and  the  powder 
when  heated  by  the  induced  RF  currents  and  cooled  by  the  walls  of  the 
container  was  completed.  We  find  that: 

°  The  molten  region  tends  to  take  the  shape  of  a  flattened  body  of 
revolution;  it  acts  as  a  fin,  with  the  source  of  heat  in  a  ring 
near  the  equatorial  plane  naar  the  surface  of  the  body.  Heat 
travels  inwards  and  is  lost  by  radiation  and  conduction  along 
the  way. 

#  If  a  substantial  heat  loss  from  the  top  surface  of  the  melt  occurs, 
the  shape  of  the  melt  is  vary  sensitive  to  the  amount  thus  lost. 

#  There  should  be  a  large  amount  of  powocr  to  insulate  the  melt  on 
its  under  side,  although  the  exact  amount  is  not  critical.  This 
means  that  the  molten  zone  should  be  confined  to  the  upper  portion 
of  the  cooled  container. 

Thermal  models  have  been  developed  for  typical  skull  melting  systems 
and  the  effects  of  hoat  sources  and  sinks  upon  melt  shape  have  been 
computer-profiled  and  are  discussed  in  Appendix  B. 

2.4  Control  of  RF  Power 


Precise  control  of  the  RF  power  is  required  for  successful  crystal 
growth.  There  is  a  basic  tendency  toward  instability  when  induction 
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heating  is  used  direcrly  to  melt  semiconductors  and  refractory  oxides. 
Therefore,  a  radio  frequency  current  control  system  will  be  used  to 
maintain  a  constant  RF  magnetic  field. 

The  balance  between  heat  loss  and  heat  generation  processes  result  in 
a  stable  sysfem  in  tnis  case.  Past  experience  shows  that  the  control 
system  response  time  can  be  as  fast  as  one  second.  The  thermei  time 
constant  of  the  charge  will  be  longer  than  this,  with  the  result  that 
temperature  fluctuations  due  to  external  disturbances  will  be  small. 

we  expect  to  obtain  fields  as  strong  as  I04  amperes  per  meter  if 
necessary,  which  will  permit  utilization  of  the  full  output  of  the 
RF  power  '-ource. 

2.5  Design  and  Construction  of  the  Cold  Crucible  Assembly 

Based  upon  the  elec+rical  and  thermal  analyses  which  have  been  carried 
out,  we  have  completed  the  design  of  a  cold  crucible  assembly  which  is 
shown  in  Figure  I.  The  copper  structure  has  been  designed  for  high 
pressure  operation  and  integra+ion  with  the  ADL  crystal  growing  furnaces 
presently  in  use  at  the  Air  Force  Cambridge  Research  Laboratories. 
Construction  of  the  cold  crucible  assembly  has  been  initiated. 
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3.0  FUTURE  WORK 


°  The  construction  of  the  cold  crucible  assembly  will  be  completed. 

0  Upon  receipt  and  installation  of  the  high  frequency  RF  power  supply 
(Heuttinger  Elektronik),  testing  and  evaluation  of  the  cold  crucible 
assembly  will  be  initiated. 

°  Performance  characteristics  of  the  cold  crucible  assembly  will  be 
determined  and  experiments  involving  the  melting  of  refractory  oxides 
under  oxidizing  conditions  will  proceed. 

#  During  the  course  of  these  experiments,  it  is  anticipated  that  the 
criteria  for  integration  of  the  cold  crucible  assembly  within  a 
controlled  atmosphere  furnace  will  be  developed. 

°  The  cold  crucible  assembly  (less  control  system)  will  be  delivered  to 
the  Air  Force  Cambridge  Research  Laboratories. 
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APPENDIX  A 


THEORETICAL  ANALYSIS  OF  THE  ELECTRICAL  CHARACTERISTICS  OF  SKULL  MELTING 


A. I  Effect  of  the  Cage 


We  assume  that  the  cage  consists  of  vertical  copper  tubes  parallel  to 
the  axis  of  the  coil  as  shown  in  cross  section  in  Figure  A. I.  The 
effect  of  the  bottom  of  the  cage  is  small  and  will  be  neglected  here. 


Coil 


Cage 


CROSS  SECTION  OF  SKULL  MELTING  CAGE 
FIGURE  A. I 


If  the  coil  ha~  n  turns  ano  a  l^ngtt  l  cr.d  carries  an  alternating 
current,  I,  the  rms  current  per  unit  length  is 

i  =  ^  (A/m)  ( I ) 

and  the  axial  magnetic  field  is 


H_  =  7—  (A/m) 
0  l 


(2) 
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At  frequencies  in  the  range  I -10  MHz,  the  skin  depth  et  penetration  of 
the  field  into  the  copper  tubes  of  the  cage  is  only  about  2  x  I0“5m. 
The  copper  therefore  essentially  excludes  the  field  and  acts  like  an 
almost  perfect  diamagnetic  material.  Under  these  conditions  each  tube 
has  a  circulating  current  per  unit  length  equal  to  -i  (A/r). 


The  energy  stored  in  the  field  is 


E-rV*  <Acon  ‘W  (Joule) 


where  the  A's  are  the  cross-sectional  areas  of  coil  and  cage  and  xnere 
Uq  =  4*  x  I0“7  Henry/m  and  is  the  permeability  of  free  space.  Fror- 
Equation  2  the  expression  becomes 


E 


%  n  1 


2'2  (A  -  A  ) 
ooil  cage 


(4) 


The  energy  can  also  be  expressed  in  terms  of  the  net  self- inductance, 
L,  of  the  assembly  by  the  relation 


E  = 


On  equating  4  and  5,  we  obtain  for  the  inductance 

(A  -  A  ) 

L  >  rA.0  0011  j  ^  (Henry) 


The  cage,  therefore,  reduces  the  inductance  of  the  coil  by  a  *azr: r 
**coil  ”  *cage),Acoi I* 

The  cage  also  adds  some  series  resistance  to  the  coii.  The  circumfer¬ 
ential  current,  i,  in  a  tube  of  the  cage  is  distribution  in  dep+h,  x, 
into  the  copper  according  to  the  approximate  formula 
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(7) 


J 


1  -x/6 

6  ® 


where  J  is  the  current  density  (A/m  ),  and  6  is  the  skin  depth,  given 
by  the  relation 


OWqTO 

Here  f  is  the  frequency  (Hz)  and  a  is  the  conductivity  of  copper 
(mhos/m) . 


The  power  dissipated  per  unit  area  of  tube  surface  is  then 


Jl 

2o6 


2 

(W/m  ) 


(9) 


The  total  power  for  N  tubes  of  height,  h,  and  diameter,  D,  is 


_  N  *0  h  i2 
2o6 


n2  Nt  Dh  |2 
2ol26 


(10) 


and  the  effective  series  resistance  introduced  into  the  coil  is 
therefore 


R  =  slmsi  (M 

2olz6 

We  r-we  assumed  that  h<t.  If  the  tubes  are  longer  than  the  coil,  h, 
in  Equation  II,  is  replaced  by  l . 

A. 2  Heat  Input  to  the  Charge 

We  assume  for  mathematical  convenience  that  the  charge  has  the  shape 
of  an  oblate  spheroid,  which  is  a  good  approximation  to  the  actual 
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disk-like  shape.  The  equation  of  the  surface  is 


pZ  ft, 

—  +  —  =  I 


where  a  and  c  are  the  semi -axes  of  the  spheroid,  and  2  and  r  are 
vertical  and  horizontal  cylindrical  coordinates. 

The  distribution  of  induced  current  in  the  spheroid  depends  strongly 
on  the  skin  depth,  6,  in  the  molten  material.  The  two  limiting  cases, 
6»c  and  S«c,  are  mathematically  tractablo,  while  the  intermediate 
case,  6*vc,  is  too  difficult  for  analytical  treatment.  The  limiting 
cases,  however,  serve  to  bracket  the  intermediate  case  and  give  con¬ 
siderable  insight  into  the  general  problem. 

A.2.1  Large  Skin  Depth 

When  6»c,  the  uniform  magnetic  field,  Hq,  produced  by  the  coil, 
penetrates  tf  •  charge  with  little  distortion.  Under  these  con¬ 
ditions  the  induced  circumferential  electric  field,  EQ,  is  given 
by  the  simple  circuital  emf  equation. 

2*r  E0  =  -  ~  Ctr2  B)  *  -  wr2  n0  ju  HQ  (13) 

where  wr2  B  Is  the  instantaneous  magnetic  flux  through  a  circle 
of  radius,  r,  and  <0  is  the  angular  frequency. 

The  associated  current  density  is 


J8  *  •  Ee  *  -  i  m0  “  r  Ho 


The  rate  of  heat  generation  per  unit  volume  at  the  radius,  r, 
is 
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=5-lJol2=?V“2 


(15) 


r2  Hq2 

The  power  per  unit  area  projected  on  the  (x,y)  plane  is 
Equation  12, 

W"»  =  2z  W,,f  =  ^  ouA2  “2  Ha2  c  r2  7 1  -  — 

2  0  0  a2 

It  is  to  be  noted  that  W"  is  zero  at  both  r  =  0  and  r  * 
is  maximum  at  r  =  a/2/3  =  0.816a. 


The  total  heating  rate  for  the  whole  ellipsoid  Is 


2tr 


W  =  /  W"  2trrdr  =  ^  owQ2  u2  a4  c  HQ2 


From  Equation  8,  this  can  also  be  written,  in  terms  of 
of  a,  as 


,,  4.  “0  “  a“  c  V 

I*  «2 


As  a  numerical  example  we  take  the  following  values: 


Pq  *  4ir  x  I0‘7  H/m 
f  =  10  x  I06  Miz 
a  *  0.04m 
c  *  0.01* 

4  *  2c  =  0.02m 
Hq  *  1000  A/m 

Then  W  »  4000  watt. 


from 

(16) 

a  and 

(17) 

i  instead 

(18) 
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The  value  chosen  for  6  corresponds  to  o  *  63  mho/m  and  Is 
probably  near  the  limit  of  validity  of  the  theory.  The  value 
of  Hq  corresponds  to  that  tor  a  5-turn  col  I  of  length  10cm 
carrying  a  current  of  20  amperes. 

A.2.2  Small  Skin  Depth 

The  case,  6«c,  may  not  occur  In  practice  for  the  range  of  con¬ 
ductivities  likely  to  be  formed  for  molten  refractory  materials, 
without  the  use  of  frequencies  well  in  excess  of  10  MHz.  How¬ 
ever,  It  is  still  of  value  to  determine  the  power  distribution 
in  the  charge  for  this  limiting  case. 

Like  the  cage  tubing,  the  spheroid  acts  like  a  nearly  perfect 
diamagnetic  body  when  d«c.  The  current  distribution  on  the 
surface  of  the  spheroid  is  such  that  the  external  magnetic  field 
is  the  same  as  that  produced  by  a  uniformly  magnetized  body  of 
the  same  shape.  The  Intensity  of  magnetization,  m,  (dipole 
moment  per  unit  volume)  is  related  to  the  applied  field,  H^, 
of  the  coi I  by  the  formula 


where  N  is  the  demagnetization  factor,  which  depends  on  the 
ellipsoid  shape.  For  a  sphere  N  *  1/3.  For  an  oblate  spheroid 
with  semi-axes,  a,  a,  c,  the  factor  is  given  by^^ 

N  *  (-——■)  (I  -  i  tan’*  k)  (20) 

k*  k 


where 
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(21) 


Table  A. I  shows  how  N  varies  with  c/a. 


TABLE  A. I 


DEMAGNETIZATION  FACTOR  N 


c 

N 

a 

0.0 

1.000 

0.2 

0.750 

0.4 

0.588 

0.6 

0.478 

0.8 

0.396 

1.0 

0.333 

For  the  shape  of  Interest  here,  where  c/e  -  1/4,  N  »  0.7037. 
For  this  value  the  intensity  of  magnetization  is  from 
Equation  19 


m  »  -  3.375  Hq  (22) 

Knowing  m,  and  therefore  the  effective  distribution  of  magnetic 
charge  density  over  the  surface  of  the  spheroid,  one  can  cal¬ 
culate  the  magnetic  field  distribution  and  from  it  the  equivalent 
surface  current  distribution. 

In  Flg-.re  A.2  we  consider  the  ellipsoid  as  a  perfect  diamagnetic 
body  with  permeability  y  *  0.  Then  the  magnetic  induction,  Bj, 
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THE  ELLIPSOID  AS  A  PERFECT  DIAMAGNETIC  BODY 


FIGURE  A. 2 


within  the  ellipsoid  is  given  by 

Bj  «  uHj  *  0  (23) 


Also, 


Bj  «  Hj  +m  (24) 

The  directions  of  Hj  and  m  are  shown  In  Figure  A.2. 

We  now  apply  the  boundary  condition  that  the  tangential  com* 

ponent  of  H  is  continuous.  Now  the  field,  H  just  outside 

the  ellipsoid  is  known  to  be  tangential,  as  indicated  in 

Figure  A.2.  Therefore,  just  below  the  surface,  Hj  »  H^,  as 

shown  In  the  vector  diagram  In  Figure  A. 3  where  the  vector 

sum  of  Ha  and  H  equals  the  Internal  field,  H. , 
on  * 


MAGNETIC  FIELD  VECTORS 


FIGURE  A. 3 


From  the  diagram  we  obtain  the  result 
He  *  Hj  sin  *  *  -  m  sin  t 
where  the  angle,  ♦,  is  defined  In  Figure  A. 2. 

On  substitution  for  m  from  Equation  15  we  find  that 

H  sin  * 
h  s  -y- 
ne  I  -  N 

The  surface  current,  (A/m),  is  therefore 

U 

,  _  0  c  cos  ♦ _ 

I  -  N  - - 

/a2  sin2  ♦  +  c2  cos2  + 


(25) 


(26) 


(27) 
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The  total  power  Is 


W  =  /  Rs  I2  2*  rds 


(28) 


where  R  Is  the  surface  resistance,  (ohm/0),  given  by  Equation 
4  as 


*sS2k 


(29) 


and  ds  is  an  element  of  arc  of  the  ellipsoid  given  by 


ds  *  /a*  sin*  4  +  c-  cos2  ♦  *  d# 


(30) 


From  Equations  27-30  we  obtain  for  the  power 
v/2 


W  = 


cos3  4d+ 


H02  j  _ 

o6(|  -  N)2  w  /a14  cos'4  ♦  +  c2  sin2  V 
'  2 


a2H  2  wc  F(f) 
o  c 

(|  -N)2 


(31) 


where 


(2x2  -  i)  log.  (/x*  -  I  +  x) ^ — 

-■^TTTiTi  • 


(32) 


Table  A.2  shows  how  the  second  factor  in  Equation  31  depends 
on  a/c.  It  can  be  seen  that  the  power  is  quite  constant  for 
c/a  between  0.1  and  1.0  (for  fixed  a).  For  very  flat  ellipsoids 
the  power  increases  slowly. 


25 


TABLE  A. 2 


GEOMETRICAL  FACTOR  IN  EQUATION  31 


C  ire  F  - 

—  c 


a 

a(l  -  N>2 

0.0001 

24.020 

0.001 

18.101 

0.01 

12.534 

0.1 

8.190 

0.2 

7.510 

0.25 

7.465 

0.4 

7.564 

0.6 

8.121 

0.8 

8.744 

1.0 

9.425 

The  radial  distribution  of  the  power  Is  also  of  Interest.  From 
Equation  28  we  find,  on  counting  the  power  on  both  top  and  bottom 
of  the  ellipsoid,  that 


dW 

dr 


4»R  H  2C2 
s  o 

(I  -  N)2 


where 


G  <-) 


(33) 


(34) 


Table  A. 3  shows  how  G  varies  with  r/a  for  the  case,  c/a  *  1/4. 
The  last  column  In  the  table  gives  the  fraction  of  the  power 
within  a  radius,  r.  It  Is  seen  that  the  power  Is  strongly 
concentrated  at  the  equator  of  the  ellipsoid. 
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TABLE  A. 3 


RADIAL  DISTRIBUTION  OF  POWER  FOR  c/a  =  1/4 

£ 

a 

G(-) 

9 

W(r) 

W(a> 

0.0 

.0000 

.0000 

0.1 

.0010 

0.2 

.0083 

.0001 

0.3 

.0296 

0.4 

.0757 

.0021 

0.5 

.1650 

0.6 

.3317 

.0121 

0.7 

.6532 

0.8 

1.3492 

.0549 

0.9 

3.4094 

1.0 

00 

1.0000 

For  c/a  =  1/4  the  total  power,  from  Equation  31  and  Table  A. 2, 
is 


W  = 


7.47 


(35) 


We  expect  that  this  expression  will  remain  approximately  valid 
for  6  *  c/2  -  0.005m.  At  f  *  10  Wz,  this  corresponds  to 
(i  3  1000  mho/m.  For  a  *  0.04m  and  Hg  =  1000  A/m,  as  before, 
we  find  from  Equation  35  that  W  =  2400  watt  which  is  to  be 
compared  with  the  value  of  4000  watt  calculated  on  the  6»c 
assumption  for  6  =  0.02m. 


On  substituting  for  o  from  Equation  8  we  find  that  Equation  35 
takes  the  form 


(36) 


W  =  3.73  yQ  u  a2  HQ2  6 

which  shows  that  W  tends  to  zero  for  very  small  6.  The  limiting 
expressions,  Equations  36  and  18,  therefore  show  that,  as  6 
varies  from  small  to  large  values,  due  to  variation  of  the  con¬ 
ductivity,  W  passes  through  a  maximum  which,  for  the  numerical 
values  chosen  above,  probably  lies  between  2000  and  4000  watts 
and  occurs  for  6  =  c  =  0.0 Im. 

A. 3  Coil  Resistance  Due  to  Load 


On  substituting  for  from  Equation  2  Into  Equation  36,  we  can  express 
W  in  terms  of  the  current  I  in  the  coll  by  the  formula 


4.78  u_  tit  a2  6  n2  I2 


W  = 


(37) 


The  load  therefore  acts  like  a  resistance 


4.78  Pq  u  a2  $  n2 


(6«c> 


(38) 


Likewise,  from  Equation  18  we  obtain 


y0  »  a-  c  n* 

Rl  ■  TT - - (s>>cl 


(39) 


A. 4  Startup 

In  order  to  start  the  inductive  heating  one  must  usually  add  some 
high-conductivity  material  such  as  particles  of  the  same  metal  as  is 
present  in  the  composition  of  the  charge  itself.  We  consider  a  metal 
sphere  of  radius  a_  and  conductivity  o.  In  the  high  frequency  field, 
Hq,  the  sphere  acts  as  a  diamagnetic  body  with  a  demagnetization 
factor  N  *  1/3  (see  Table  A. I).  From  Equation  19  the  effective 
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Intensity  of  the  magnetization  of  the  sphere  Is 


m  =  -  j  HQ  (40) 

The  total  magnetic  moment  of  the  sphere  Is 

M  =  ^wa3m  =  -  2»a3  Hq  (41) 

For  the  special  case  cf  the  sphere,  unlike  that  of  a  spheroid,  the 
moment,  M,  acts  as  though  It  were  located  at  the  center  of  t.ie  sphere. 
The  tangential  magnetic  fiald  at  the  surface  of  the  sphere,  due  to  M, 
is 


Cho‘dlpol.  •  ■ ^ ■  - *  *  s,n  6 

where  6  is  the  polar  angle. 

Since  the  applied  field,  Hq,  gives  a  tangential  component,  ~  HQ  sin  e, 
at  the  surface  of  the  sphere,  the  total  tangential  field  is 

H0  *  "  I  *\)  sln  6  (43) 

The  circumferential  current  per  unit  length  on  the  sphere  is  therefore 

*  *  I  H0  s,n  0  (44) 

since  it  equals  the  discontinuity  in  (!q  at  the  surface  of  the  sphere. 

The  surface  resistance  is  given  by  Equation  32.  On  multiplying 
this  by  i^2  and  integrating  over  the  area  of  the  sphere,  we  obtain 
for  the  power  input  to  the  sphere 
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It 

W  -  / 
0 


(2^6H!H02  Sl°2  0)  2l,a2  Sln  6d8 


Ho2a2 

06 


(45) 


It  Is  of  Interest  to  consider  that  this  power  is  balanced  by  blackbody 
radiation  from  the  surface  of  the  sphere  at  absolute  temperature,  T, 
with  the  surroundings  at  temperature,  T.  Then 

3irH  2 

W  o0(T»  -  T0*>  -  -^2-  (46) 


or 


°o<T4  - 


3Ho2 


(47) 


where  cQ  =  5.672  x  I0"8  W  m"2  deg-1*  is  the  Stefan-Bol tzmann  constant. 
Equation  47  shows  that  the  equilibrium  temperature  is  independent  of 
the  radius  of  the  sphere.  For  the  case  of  a  field,  HQ  =  1000  A/m, 
as  before,  and  with  a  =  5.0  x  I0“5m,  oQ(T4  -  TQ4)  3  1500  W/m  ,  and 
T  -  431  K.  The  metal  sphere  would  therefore  not  melt  for  the  field 
assumed.  However,  for  a  ten  times  greater  field  of  10,000  A/m,  which 
is  not  unreasonable,  we  get 


and 


OqCT4  -  y)  *  150,000  W/m2 

T  -  1276  K 


which  will  melt  certain  metals. 

A  more  appropriate  assumption  for  skull  melting  is  that  the  metal  sphere 
is  embedded  in  a  powder  of  thermal  conductivity  K.  Then  the  equilibrium 
temperature  is  given  by 
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(48) 


T  _  T  =  _ _ 

'  'o  4*Ka 

For  K  =  0.1  watt/(m  deg  K),  a  =  !  x  I0"3m,  Hq  =  10,000  A/m,  and  tne 
other  constants  the  same  as  before,  we  find  that  M  *  1.88  watt  and 
T  =  I600°K.  The  particle  would  therefore  raise  the  temperature  of 
the  surrounding  powder  to  a  value  where  Inductive  heating  of  the 
powder  would  occur. 
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APPENDIX  B 


THEORY  OF  THE  THERMAL  CHARACTERISTICS  OF  SKULL  MELTING 

It  has  been  established  already  that  the  currents  produced  In  an  oblate 
spheroid  by  an  induction  heater  are  largely  concentrated  around  an  equa¬ 
torial  ring  which  is  located  near  the  surface.  While  the  current  density 
as  a  function  of  position  has  been  derived  rigorously  only  for  the  oblate 
spheroid,  the  results  for  this  case  give  us  much  insight  and  enable  us  to 
anticipate  what  will  probably  happen  in  the  case  of  a  more  complicated 
body  of  revolution  whose  radius  is  large  compared  with  its  thickness  in 
the  axial  direction  and  which  is  an  electrical  conductor. 

Qualitatively  speaking,  we  may  say  that  the  currents  in  the  body  will  estab¬ 
lish  themselves  in  such  a  way  as  to  exclude  most  of  the  field  impressed  by 
the  induction  heater  coil.  The  value  of  B  produced  by  a  ring  of  current 
depends  mainly  on  the  current,  and  not  so  much  on  the  radius  of  the  ring 
(in  a  solenoid,  the  value  of  H  is  independent  of  the  radius  of  the  coil), 
whereas  the  total  flux  produced  is  the  product  of  the  area  of  the  ring  and 
the  average  value  of  B  within  the  ring.  Therefore,  in  the  case  of  a  con¬ 
ductor  which  is  a  flattened  convex  body  of  revolution,  one  may  expect  that 
the  induced  current  will  be  largely  concentrated  in  a  ring  of  the  largest 
possible  radius  which  Is  contained  in  the  body.  It  is  probably  best  to 
leave  this  kind  of  qualitative  argument  in  its  present  rather  vague  state, 
and  simply  to  conclude  that  It  is  of  Interest  to  study  the  shape  of  the 
isotherm  at  the  melting  temperature  In  a  region  filled  with  a  powder  of 
low  thermal  conductivity,  which  contains  a  ring  source  and  which  is  fairly 
closely  contained  by  a  vessel  (represented  by  an  isotherm  at  300°K)  main¬ 
tained  at  room  temperature.  We  will  also  suppose  that  near  +he  top  of  the 
vessel  a  rather  large  amount  of  heat  Is  lost  in  order  to  simulate  the  con¬ 
ditions  that  must  obtain  when  the  crystal  is  being  pulled  slowly  from  the 
melt  through  an  aperture  above. 


32 


One  way  to  approach  the  problem  Is  to  take  the  shape  of  the  300K  Isotherm 
from  the  drawing  of  the  tubes  containing  the  coolant  and  to  solve  rigorously 
the  boundary  value  problem  thus  defined.  However,  the  physical  parameters 
involved  (for  example  the  thermal  conductivity  of  the  powder,  which  will  de¬ 
pend  upon  its  packing  and  any  sJntering  that  may  take  place)  are  very  poorly 
known,  so  it  is  not  justified  to  use  such  rigorous  procedure.  Instead  we 
may  employ  a  simpler  technique  and  obtain  at  least  a  preliminary  understanding 
of  the  effects  on  the  molten  zone  that  will  be  produced  by  variations  in  the 
main  features  of  the  system.  Without  great  difficulty  some  important  results 
can  thus  be  obtained  which  will  serve  as  guide  posts  in  the  experimental 
program. 

In  doing  the  simplified  problem  we  are  at  liberty  to  place  sources  and  sinks 
of  various  strengths  at  such  locations  as  we  choose,  provided  that  these 
locations  are  outside  the  region  of  interest,  which  we  will  take  to  be  bounded 
by  the  3Q0K  isotherm.  Of  course  the  ring  source,  which  is  used  to  model  the 
induction  heating,  must  be  within  the  melt.  We  make  use  of  external  sources 
and  sinks  to  model  the  loss  of  heat  from  the  top  and  to  control  the  shape  of 
the  300K  isotherm. 

It  appears  that  the  3Q0K  isotherm  must  pass  close  to  the  ring  source  in  the 
melt  in  order  to  gUe  the  isotherm  at  the  melting  point  its  expected  flat 
shape.  Also,  we  want  the  300K  isotherm  to  be  close  to  the  melt  at  the  top  (to 
model  the  loss  of  heat  through  the  aperture  used  to  extract  the  crystal)  and 
otherwise  to  resemble,  at  least  roughly,  the  shape  of  the  cup  formed  by  the 
tubes  containing  the  coolant. 

In  Figure  B. I  we  show  what  we  have  taken  to  be  the  "standard"  arrangement  for 
the  300K  and  2318K  isotherms  (the  melting  point  of  alumina  is  23I8K).  The 
positions  of  the  sources  and  three  sinks  are  shown  and  the  drawing  is  to  scale, 
as  indicated.  The  arrangement  pictured  in  Figure  B. I  was  obtained  after  a 
number  of  mathematical  experiments;  one  lesson  that  was  learned  very  early 
is  that  the  300K  sink  on  the  axis  has  a  drastic  effect  on  the  shape  of  the 
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melt,  and  that  the  lower  surface  of  the  melt  needs  to  be  well  Insulated. 

The  melt  Is,  In  effect,  acting  as  a  fin,  with  the  heat  traveling  towards  the 
center  nnd  being  lost  alonq  the  way.  The  Insulation  along  the  under  side, 
which  serves  to  keep  the  heat  loss  to  a  minimum,  can  be  obtained  by  making 
sure  that  the  molten  zone  Is  In  the  top  part  of  the  cooled  container,  so  that 
a  thick  blanket  of  powder  is  left  unmelted  beneath  the  231 8K  Isotherm. 

We  will  now  illustrate  the  sensitivity  of  the  shape  of  the  molten  region  to 
changes  In  some  of  the  system  parameters  by  means  of  figures  which  show  the 
two  standard  Isotherms  and  two  slightly  distorted  Isotherms.  The  distorted 
isotherms  are  obtained  by  solving  the  problem  with  some  changes  in  the 
boundary  conditions. 

In  Figure  B.2  we  show  two  pairs  of  Isotherms  to  demonstrate  the  sensitivity  of 
the  shape  of  the  melt  to  a  variation  in  the  position  of  the  upper  part  of  the 
3Q0K  isotherm  from  the  standard  position.  The  lower  part  of  the  300K  Isotherm 
has  been  held  nearly  in  place.  It  is  clearly  seen  that  the  top  surface  of  the 
melt  moves  further  than  the  300K  isotherm,  while  even  the  bottom  part  suffers 
a  significant  displacement.  It  is  of  interest  to  "ote  that  the  changes  seen 
are  due  mainly  to  a  change  of  50  watts  in  the  strength  of  the  sink  on  the  axis 
(from  300W  to  250W). 

In  Figure  B.3  we  show  two  pairs  of  closely  spaced  Isotherms.  The  lower  part 
of  the  300K  isotherm  has  been  caused  to  move  about  3mm  below  its  standard 
position  (on  the  axis),  but  the  corresponding  variation  in  the  position  of  the 
23I8K  isotherm  is  only  about  Imm  at  the  bottom  and  nearly  nothing  at  the  top 
of  the  melt.  Thus  it  Is  seen  that  the  shape  of  the  molten  zone  is  not  very 
sensitive  to  variations  in  the  geometry  of  th»  uottom  of  the  cup,  provided 
that  the  blanket  of  insulation  is  present. 

We  also  determined  the  sensitivity  of  the  shape  of  the  molten  zone  to  a  change 
in  the  input  power  (supplied  by  the  induction  heater)  by  raising  the  strength 
of  the  ring  source  from  3150  W  to  3300W,  while  the  shape  of  the  300K  isotherm 
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was  kept  essentially  the  same  as  the  standard.  The  result  is  shown  in  Figure 
B.4,  which  indicates  that  the  shape  of  the  melt  is  only  slightly  affected  by 
the  change  in  power  input.  This  is  probably  due  to  the  fact  that  the  distance 
between  the  side  of  the  molten  zone  and  the  cold  container  is  very  small,  so 
a  minute  change  in  this  dimension  is  all  that  is  needed  to  cause  most  of  the 
added  power  to  flow  directly  into  the  container  walls  near  the  plane  z  =  0. 

Under  certain  circumstances  it  might  be  desirable  to  let  the  molten  mass  cool 
slowly,  of  its  own  accord,  without  penetrating  the  layer  of  powder  that  lies 
over  the  melt.  In  such  a  case  one  might  obtain  a  symmetrical  melt,  as  shown 
in  Figure  B.5. 

So  far,  we  have  discussed  only  the  case  of  molten  alumina.  Another  material 
of  interest  is  zirconia,  whose  melting  point  is  2988K.  In  Figure  B.6  we  show 
a  set  of  sources  and  sinks  and  the  corresponding  isotherms  at  2988K  ard  300K. 
Because  of  its  higher  temperature  the  upper  surface  of  the  pool  of  molten 
zirconia  may  be  expected  to  lose  much  more  heat  by  radiation  than  is  the  case 
with  alumina.  To  allow  for  this,  the  strength  of  the  point  sink  has  been 
set  at  800W,  anc  the  strengths  of  tne  other  sources  and  sinks  adjusted  in  an 
attempt  to  produce  a  suitable  configuration  of  the  2988K  and  300K  isotherms. 
While  the  amount  of  mathematical  experimentation  with  the  model  of  the  higher 
melting  point  material  has  been  very  limited,  we  expect  that  it  would  be 
necessary  to  separate  the  ring  source  slightly  further  from  the  main  ring 
sink  in  order  to  produce  a  2988K  isotherm  that  is  similar  in  appearance  to 
the  231 8K  isotherm  of  Figure  B. 1 .  From  an  operational  point  of  view,  this 
means  that  the  frequency  of  the  induction  locator  must  be  set  so  as  to  obtain 
a  large  skin  depth  in  the  melt. 

We  expect  that  the  sensitivity  of  the  shape  of  the  zirconia  melt  to  changes 
in  the  external  configuration  will  be  qualitatively  the  same  as  in  the  case 
of  the  alumina  melt  already  discussed. 
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FIGURE  But  ISOTHERMAL  PROJECTION 
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FIGURE  B J3  ISOTHERMAL  PROJECTION 
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APPENDIX  C 


TEMPERATURE  DISTRIBUTION  DUE  TO  A  RING  SOURCE 


In  this  section  we  derive  the  formula  for  the  temporal ure  (In  steady  state) 
at  a  point  In  an  Infinite  medium  which  contains  a  source  of  heat  In  the  form 
of  a  ring.  We  assume  the  thermal  conductivity  of  the  medium  to  have  the 
fixed  value  K  (watts  cm"*  deg  C"*),  so  the  temperature  must  satisfy  Laplace's 
equation 

<t2T  =  o 

and  furthermore  we  require  that  the  temperature  be  zero  at  Infinity. 


In  accordance  with  these  assumptions  we  find  the  temperature,  dT,  at  a  field 
point,  P,  whose  coordinates  are  (x,  y,  z)  produced  by  a  point  source  of 
strength  dQ  watts  locc+ed  at  P'  with  coordinates  (x',  y',  z')  to  be 


where 
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(I) 


s  =  /(x  -  x')*  +  (y  -  y'J^  +  (z  -  z*)^ 


(2) 


By  adding  up  the  contributions  due  to  many  such  sources  we  can,  at  least  in 
principle,  find  the  temperature  due  to  any  distribution  of  heat  sources. 


In  Figure  C. I  we  show  a  ring  souroe  of  radius,  a,  located  in  the  z-plane  and 
a  field  point  at  an  arbitrary  point.  It  is  evident  from  the  diagram  that 

x'  ■  a  cos  6 
y'  =  a  sin  0 

z'  =  0  (3) 

Now  suppose  the  souroe  to  dissipate  Q'  watts  cm"*  uniformly  along  its  cir¬ 
cumference.  Then 
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dO  =  O'  ad6  (4) 

since  ad6  is  the  element  of  arc.  Therefore  we  find  from  the  above  equation 


dT 


_ Q*  a  d8 _ 

4irK/(x  -  a  cos  8)a  +  (y  -  a  sin  0)2  +  zz 


(5) 


or,  by  integration  of  Equation  5  around  the  loop, 

T  =  ?  —  -  d9  - —■■■■■■  (6) 

*  0  Ax  -  a  cos  8)z  +  (y  -  a  sin  e)2  +  z2 

On  account  of  the  symmetry  of  the  problem  we  can  replace  x  in  Equation  6  by 
r  and  replace  y  by  0.  Thus 

T  .  Sll  f . . __g . . --  O) 

0  i/r2  +  a2  +  z2  -  2ar  cos  8 


The  integral  in  the  right-hand  side  of  Equation  7  is  not  "elementary",  but 
can  be  evaluated  simply  in  terms  of  the  well-known  complete  elliptic  Integral, 
Ktm),  which  Is  defined  by  the  formula 

w/2  ha 

K(m)  =*  /  qT  ■  ■■  —  C8) 

0  /I  -  m  sin2  4 


First  make  the  substitution 


8  =  24  +  * 


In  Equation  7  to  obtain 


T  .  21»  ... _ 

-w/2  Aa  +  r)2  +  z2  -  4ar  sin*  ♦ 


Jit 


C9) 


Next  remove  the  factor  (a  +  r)2  +  z2  from  the  expression  whose  square  root 
is  indicated  and  note,  since  sin2  4  is  even,  that  the  range  of  integration 
can  be  reduced  by  setting  the  lower  limit  to  zero  if  a  factor  of  two  Is 
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placed  In  front  of  the  integral.  Thus  we  find 


T  = 


irK/(a  +  r)*  +  zz  0  / 


(10) 


1  ~  (a  ■*-  ^  »'*  ♦ 


Whence,  by  the  definition  of  Equation  8,  it  is  clear  that 


T  » 


Pla _ 7t  r  4ar 

—  Sa  +  r)*  +  zz*^ 


irK^a  +  r)z  +  zl 


(II) 


The  function  K(m)  Is  easily  computed  by  means  of  the  approximations  given  by 
Hastings.*  Finally  we  let  Q  be  the  total  heat  dissipated  by  the  ring  so 


Q»  =  Q/2»a 


and 


T  = 


2irzK/(a  +  r)2  +  zk 


-  t-7aT 


4ar 

~r)r 


+  2r 


(12) 


If  the  ring  is  located  at  a  height,  h,  instead  of  at  the  plane,  z  -  0,  the 
quantity,  z,  In  Equation  12  must  be  replaced  by  the  quantity,  z  -  h.  Note 
that  K(0)  =  w/2,  so  a  ring  source  of  radius,  a  *  0,  produces  a  temperature 
field 

T  - - (13) 

4wK/rz  +  zz 

This  is  obviously  the  field  produced  by  a  point  source,  as  it  jolly  well 
ought  to  be. 


*C.  Hastings,  Jr.:  Approximations  for  Digital  Computers,  Princeton  Univ¬ 
ersity  Press,  Princeton,  New  Jersey. 
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FIGURE  C.1  GEOMETRICAL  ARRANGEMENT  OF  THE  RING  SOURCE  A  FIELO  POINT 
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APPENDIX  D 


ANALYSIS  OF  RADIO  FREQUENCY  POWER  CONTROL 

The  radio  frequency  power  used  to  heat  the  charge  will  be  controlled  by  a 
feedback  system  monitoring  the  output  RF  current.  Controlling  the  current 
controls  the  RF  magnetic  field  and  leaves  the  conductivity  of  the  charge  as 
the  only  uncontrolled  variable  determining  the  power  Into  the  charge.  Con¬ 
ductivity  Increases  with  Increasing  temperature.  For  small  penetration  depths 
increased  conductivity  results  in  a  decrease  In  input  power  so  that  the  system 
is  unconditionally  stable.  For  large  penetration  depths  the  power  input  is 
proportional  to  conductivity,  which  varies  at  most,  approximately  as  the  first 
power  of  temperature.  The  power  loss  from  the  charge  is  mostly  by  radiation 
and  thus  varies  as  the  fourth  power  of  the  absolute  temperature.  Therefore, 
the  system  will  be  stable  for  large  skin  depth  also. 

A  typical  current  feedback  control  system  is  shown  in  Figure  D.l.  The  H-Azar 
recorder  is  used  as  a  preamplifier,  the  Input  summing  point,  and  also  to  re¬ 
cord  the  error  signal.  The  characteristic  time  constant  of  the  H-Azar,  Ti, 
is  0.5  second.  The  time  constants  of  the  RF  generator  are  approximately  this 
same  value. 

The  integral  (reset)  action,  kc,  of  the  series  00  controller  is  made  as  large 
as  possible  to  obtain  accurate  steady-state  control,  while  the  gains  of  the 
H-Azar  and  the  series  80  controller,  k|  and  kj,,  respectively,  and  the  deriva¬ 
tive  (rate)  action,  ka,  of  the  series  80  controller  are  adjusted  to  obtain 
stable  control  with  minimum  response  time. 

The  measured  frequency  response  of  a  typical  control  system  is  shown  in  Figure 
D.2.  The  l/f  slope  In  the  low  frequency  end  of  the  open  look  response  is  due 
to  the  integral  action  of  the  controller.  This  system  has  been  adjusted  for 
maximum  speed  of  response,  which  resulted  In  the  60°  phase  margin,  the  minimum 
acceptable  in  good  control  practice.  The  resultant  closed  loop  response  is 
good  to  about  I  Hz. 
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The  last  factor  In  the  system  dynamics  ts  the  thermal  transfer  function  of 
the  charge.  We  can  estimate  the  thermal  time  constant,  T3,  by  assuming  that 
the  heat  loss  Is  entirely  due  to  radiation  and  using  the  relationship 


The  heat  capacity,  C^,  =4.18  Vds  (joules/degree)  In  which  V,  d  and  s  are  the 
volume,  density  and  specific  heat  of  the  charge,  respectively.  Typical  values 
are 

V  =  40  cm3 
d  =  5.5  gm/cm3 
s  =  0. I  cal/gm°K 

which  result  In  Ch  =  92  joule/°K 

The  heat  loss  due  to  radiation  Is 

Wr  =  5.67  x  |0"12  S  T4  (watts) 

In  which  S  Is  the  radiating  surface  area 

dWr/dT  =  22.7  x  10  12  S  T^  (watts/degree) 

2 

Evaluating  this  derivative  for  S  =  20  cm  (approximately  one  fourth  of  the  , 
total  surface)  and  T  =  3000°K,  we  obtain  dWr/dT  =  12  watts  per  degree  and 
x3  =  8  seconds.  Since  this  is  appreciably  longer  than  the  response  time 
of  the  closed  loop,  temperature  fluctuations  will  be  small. 

One  final  consideration  is  the  maximum  available  power.  Power  is  determined 
by  the  magnetic  field  and,  therefore,  by  the  RF  current.  The  current  Is 
determined  by  the  Inductive  reactance  of  the  coll  and  the  output  voltage  of 
the  RF  generator.  Since 

H  *  r-  •  1  -  §r  a"d  L  *  "2 
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decreasing  the  number  of  turns,  n,  in  the  coil  increases  the  current,  I,  as 
n"2  and  increases  the  magnetic  field,  H,  and  consequently  the  output  power. 


Our  preliminary  coil  design  calls  for  a  four-turn  coil,  1.5  inches  long  and 
1.75  inches  mean  radius.  The  coil  inductance  is  given  by 

L  =  ToFTl?  =  10  x  1.5  4~9  x  1.75  "  1,6  rnicrohenrles 

The  factor  for  loss  of  inductance  due  to  the  cross-sectional  area  of  the 
cage,  £-9.f,  is  estimated  to  be  about  .8,  so  that  the  net 

!nductance  is  about  1.3  microhenries.  Assuming  that  the  maximum  voltage 
is  5000  volts  and  at  a  frequency  of  8  megahertz,  we  have 

.  5000 

I  =  - = - rrr - r-= - =  80  amperes 

2n  x  8  x  10°  x  I . 3  x  10  6 

and 

H  =  | . 5 »  x423?x  T0~2'  =  8000  amPeres/meter. 


Substituting  this  value  Into  Equation  17  from  Appendix  A, 


V#  *  T5  o  no2  w2  a4c  Ho2  , 

with  the  other  parameters 

a  =  50  mho/meter 
no  =  4it  x  I0"7  henry /meter 

W  =  5.0  x  I07  radians/second 
a  =  .032  meter 
c  =  .01  meter 

we  obtain 


W  =  50  Kilowatts 


which  is  the  rated  capacity  of  the  RF  generator. 
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L  &  N  Series  80 
C.A.T.  Controller 


Closed  Loop  Phase 


FIGURE  D.2  MEASURED  FREQUENCY  RESPONSE  OF  A  TYPICAL  CONTROL  SYSTEM 


